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Wo report measurements of branching fractions for the decays B Plvii, where P are the 
pseudoscalar charmless mesons 7r~, tt", 77 and 77', based on 348fb~^ of data collected with the 



4 



BABAR detector, using B and B mesons found in the recoil of a second B meson decaying as 
B D^*'^ £vi. Assuming isospin symmetry, we combine pionic brancliing fractions to obtain B{B'^ — > 
'K~ t^vi) — (1.54± 0.17(stat) ±0.09(syst)) X 10~*; we find 3.2(t evidence of the decay B'^ — > •qt'^vi and 
measure its branching fraction to be (0.64 ± 0.20(stat) ± 0.03(syst)) x 10"*, and determine B{B^ 
r^'t^ut) < 0.47x 10"'' to 90% confidence level. Using partial branching fractions for the pionic decays 
in ranges of the momentum transfer and a recent form factor calculation, we obtain the magnitude of 
the Cabibbo-Kobayashi-Maskawa matrix element \Vub\ ~ (4.0±0.5(stat) ±0.2(syst)lo'5(thcory)) x 10"'^. 

PACS numbers: 13.20.He, 12.15.Hh, 12.38.Qk, 14.40. Nd, 14.40.Aq 



The determination of the magnitude of the Cabibbo- 
Kobayashi-Maskawa matrix [l[ element \Vub\ provides a 
critical constraint on the Unitarity Triangle; the decay 
b — > uivi is a theoretically and experimentally robust 
means of measuring | V^b | ■ In the measurements described 
in this Letter, we reconstruct the b — + uii^i decay exclu- 
sively, measuring branching fractions for the processes 
S° TT-£+iye 2] and B+ ir^l+vi. These are selected 
in the recoil of the semileptonic decay B D'^*'^Ivi, 
which provides a measurement complementary to other 
studies U 01 ; this measurement is significantly more pre- 
cise than previous measurements of its kind [3, Q . Addi- 
tionally, branching fractions for the decays S"*" r]£~^vi 
and _B+ r]'£~^vi are measured, which provide potential 
additional means of determining \Vub\ as well as a probe 
into the dynamics of the jy-jy' meson system 6^ . 

We use a sample of 383 million BB pairs, correspond- 
ing to an integrated luminosity of 348 fb^^ recorded on 
the T(45') resonance by the BaBAR detector at the PEP-II 
asymmetric-energy e^e^ storage rings. The BaBAR de- 
tector provides neutral and charged particle reconstruc- 
tion and charged particle identification, and is described 
in detail elsewhere We also use a detailed Monte 
Carlo simulation (MC) fsj to estimate signal efficiency 
and signal and background distributions. 

We tag B mesons decaying as i? — > D'-*^£i'£ through the 
full hadronic reconstruction of and mesons; D*^ 
mesons are reconstructed through K~Tr~^, K~TT~^7T^Tr^ , 
K^tt^tt'-' and K'^7t^7t~ decays, and mesons through 
if ^7r+7r+ and K'^tt'^ decays; candidates are recon- 
structed as ifg — > 7r"'"7r~; and neutral pions are recon- 
structed as TT^ — *■ 77 with the requirement 115 < m^^ < 
150MeV/c^. Masses of D candidates are required to be 
within 2.3(7 of their nominal value, where the mass resolu- 
tion a ranges between 5.7 and 19.1 MeV/c^, depending on 
the decay channel; we also use a "sideband" sample of D 
candidates with reconstructed mass in a range (typically 
4cr to 7a) off the appropriate nominal mass. We require 
charged daughters of the D candidate to originate from 
a common vertex. We reconstruct D*'^ mesons as D^tt"*" 
and D+TT*^ and mesons as D°tt° and D^'y. The mass 
difference between the D* candidate and its D daughter 
must be within 3.7(7 of its nominal value; the resolution 
(7 of this difference ranges between 0.9 and 5.7MeV/c^, 
depending on the decay mode. 

Candidate mesons are paired with tracks identi- 



fied as leptons with absolute momentum \pi\ > 0.8 GeV/c 
[9]. If a 13 candidate (its daughter kaon) is charged, it 
is required have charge opposite to (same as) that of the 
corresponding lepton. The Y = D*£ system is required 
to have invariant mass my > 3 GeV/c^ and originate 
from a common vertex. Photons consistent with origi- 
nating from bremsstrahlung from this lepton or the de- 
cay D^*") £'7(7) are added to the Y system. Assuming 
that the B Yv decay hypothesis is correct, the angle 
9 BY between the directions of the (measured) Y and its 
parent B is described by 



cos Uby 



2EbEy 



2|p*B||py 



(1) 



where Eb, jtib and \pb\ {Ey , my and |pV|) are the en- 
ergy, mass and absolute momentum of the B meson [Y 
system); for the B meson, these are inferred from ini- 
tial beam energies. If the B — ^ Yv hypothesis is cor- 
rect, we have |cos0By| < 1 up to resolution; because 
cos 9by is strongly correlated with our discriminating 
variable cos^ (jjB , we impose the loose requirement that 
I cos^sy I < 5. 

To suppress background from non-BB events, we reject 
events for which the ratio of the second and zeroth Fox- 
Wolfram moments fl^ is greater than 0.5. We also reject 
events containing lepton pairs kinematically and geomet- 
rically consistent with having originated from the decay 
of a J/t/j meson. We reject D^*'>£ candidates for which 
the event contains any K'^ — > tt+tt^ candidates not over- 
lapping this D'^*H system. We require exactly one addi- 
tional lepton with absolute momentum \pi\ > 0.8GeV/c 
in the event. If the two leptons are an e"'"e~ pair, we 
require them not to be consistent with originating from 
7 — > e^e^ conversion. This second lepton is paired with 
remaining tracks (assumed to be pions), neutral pions 
and photons in the event to form B Pf-ve candidates, 
where P is one of the mesons tt^, tt*^, rj or rj'. For 
B TT^£i'g candidates, the lepton and pion are required 
to have opposite charge. B — > ■7T°£ve candidates are sub- 
ject to the additional requirement |pVo > 2.6 GeV/c, 
where jp^o | is the absolute momentum of this tt" can- 
didate. For B ~f ri£h'i candidates, 77 mesons are recon- 
structed through decays to 77, 7r"'"7r~7r'^ and tt'^tt^tt'^, with 
invariant mass requirements 500 < rrij-y < 570, 530 < 
"iTTTTir < 560MeV/c^. Charged pions from 77 tt+tt^tt^ 
decays are required to come from a common vertex; the 



5 



TT candidates are required to have absolute laboratory 
frame momentum greater than 280MeV/c (180MeV/c) 
when coming from vt+tt^tt'' (tt^tt^tt'') candidates. The 
rj' meson in B —^ rj' ivi decays is reconstructed through 
its decay 77' — > yyTr+Tr" with the 77 candidate selected as 
above; the additional pions are required to originate from 
a common vertex, and the r^Tr+Tr" system is required to 
have invariant mass between 920 and 970MeV/c^. For 
decays [P ~ tt",??, t]'), the leptons in an event are 
required to have opposite charge. 

We define the X as a charmless meson tt^, tt", 77 or 
•q' and corresponding lepton (including photons consis- 
tent with having originated from bremsstrahlung from 
it); (^BX is defined analogously to Oby] we require 
I cos 6*5x1 < 5. For each D'^*h-Pt candidate, we re- 
quire that there be no additional tracks in the event and, 
for hypothesized B°B° {B+B") events, at most 140 MeV 
(70MeV) of neutral energy (i.e., photon candidates) not 
associated with the D^*^^ or PI candidates. In the case 
that more than one pair fulfills all requirements 

for a given event and P mode, the candidate is chosen by 
smallest |cos6'_By|, then by largest absolute P momen- 
tum. Signal events containing accepted D^*H-Pi candi- 
dates have, on average, between 1.15 and 1.39 of them, 
depending on P. 

Signal yield is extracted independently for each P; 
while we implicitly allow an event to be reconstructed in 
multiple P modes, we find the induced pairwise statistical 
correlations between our measured branching fractions to 
be negligible. The signal yield is extracted through the 
quantity cos^ (pB, where is the angle between the di- 
rection of either B and the plane containing the X and 
Y momenta: 

2 , cos'^ 9by + "^cosjcosOby cosObx + cos^ 9bx 

cos (pB = n 

sin 7 

(2) 

where 7 is the angle between the X and Y momenta. For 
correctly reconstructed signal events, we have cos^ (t>B < 
1 up to resolution. 

For a i? — > P£ve decay, is defined as the squared 
invariant mass of the lepton-neutrino system, and is cal- 
culated in the approximation that the B is at rest, i.e., 

= {rriB — EpY — Ippp, where Ep and pp are, respec- 
tively, the energy and momentum of the P meson. The 
data are divided into three bins: < 8, 8 < < 16 
and > 16GeV^/c^, in each of which the yield is ex- 
tracted separately, except in the B+ -q't^v^ mode, in 
which, due to a lower reconstruction efficiency, the yield 
is measured in a < 16GeV^/c^ bin and over the full 
q^ range. The data is described as a sum of three contri- 
butions, dN/dcos^ (j)B = Nsie'Psig + NbgPhg + NcmbVcmh, 
where these Ni and Vi are the yield and probability den- 
sity functions (PDF) of: signal ("sig"), background with 
correctly reconstructed D"'^ mesons ("bg") and back- 
grounds with combinatoric D'^'^ candidates ("cmb"). 




cos^ 



°°140P - 




- W f 
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cos^ 



FIG. 1: Distributions of cos^ (I>b for B° ^ tt-I+vi (a), B+ 
TT°£+ui (b), B+ ril+ue (c) and B+ -qtvi (d) candidates; 
filled and hollow circles represent D mass peak and sideband 
data, respectively. The curves are stacked fit results for "cmb" 
(dotted), "bg" (dashed) and "sig" (solid) PDFs, as defined in 
the text. The fits are performed in bins of but are here 
shown in the full range. 



The signal PDF, Psig, is modeled as a threshold func- 
tion (constant between zero and unity, vanishing else- 
where) with finite resolution and an exponential tail (four 
parameters). The correct D background PDF, T'bg, is 
modeled as an exponential with a nonnegative constant 
term (two parameters); the combinatoric D background, 
'Pcmb, is modeled by a second order polynomial (two pa- 
rameters). These eight PDF shape parameters and the 
Vi are determined via simultaneous unbinned maximum 
likelihood fit (see Figure [T]) of dN /dcos^ (f>B to the data, 
T'sig to MC signal events, Vhg to MC background events 
(with correctly identified D'^'^ mesons) and Vcmh to the 
sideband sample. The combinatoric yield A^cmb is further 
constrained, up to statistical accuracy, by the number of 
events in the sideband sample. Total signal yields are 
found to be 150 ± 22, 134 ± 20, 55 ± 15 and 0.6 ± 3.9 
events for n^£iy£, n^ii^i, rjlivi and rj' £vi respectively. 

The B D^*'>ii'£ reconstruction efficiency is deter- 
mined via an analogous cos^ (j)B study on "double tag" 
events, i.e., events reconstructed as BB with both B 
mesons decaying as i? ^ D'^*^Ivi. The B — > Pivt re- 
construction efficiency for each q^ bin is determined from 
the MC signal sample, as are bin-to-bin migrations due to 
the finite q^ resolution, which are small (< 9%). Overall 
efficiencies, including branching fractions and reconstruc- 
tion efhciency of the recoil i?, are found, in units of 10"'^, 
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to be 1.4, 1.8, 1.1 and 0.22 for B Tv^ivi, B tt^Ivi, 
B rjii'e and B — > 77'^^^ respectively. 

Systematic uncertainties associated with physics mod- 
eling are evaluated by determining the change in the mea- 
sured branching fraction after varying independently in 
MC within current knowledge: B —^ {p, uj^iv^ branching 
fractions, B —* ir^'^tvi branching fractions, B — > 
branching fractions, the total B charmless semileptonic 
decay branching fraction, the B charmless semileptonic 
decay spectrum [iJ] , B — > Piv^ form factors (comparing 
the model by Ball and Zwicky [l2| to that of Scora and 
Isgur [l3|) and several B — > Divg branching fractions; 
the largest is found to have an effect four times smaller 
than the statistical uncertainty. We also apply uncer- 
tainties derived from those on 77 and 77' decay branching 
fractions. 

We estimate the systematic uncertainty associated 
with the accuracy of BB background simulation by com- 
paring the cos^ 4>B distributions in signal-depleted data 
and MC samples. From study of 37 fb"""^ of e+e~ colli- 
sions 40MeV below the T{AS) resonance, we determine 
that there is no contribution from non-SB events to the 
signal; the precision to which this can be determined is 
also taken as a systemic uncertainty. 

Final state radiation in -K^t^vi decays is deter- 

mined, from simulation, to cause bin migrations no 
greater than 1.2%, which is conservatively applied as a 
systematic uncertainty, as well as to the other branching 
fractions. We apply a 0.59% (1.7%) systematic uncer- 
tainty for B^B^ {B^B~) decays associated with the as- 
sumption that double tag events can be used to estimate 
the single tag efficiency reliably. 

As double tag events are used to determine the D^*'>£:yi 
reconstruction efficiency, detector simulation uncertain- 
ties are apphed only to particles on the P£ side: 0.36% 
per track, 3% per 7r°, 2% (3%) per electron (muon). 
There is a 1.1% systematic uncertainty from counting 
BB pairs and a 1.4% systematic uncertainty from 
the T(4S') B°B° fraction Measured branching 

fractions and associated uncertainties are given in Table 
m Quoted statistical uncertainties are due to the finite 
size of data and MC samples. We combine l^vg 
and B~^ — > ir^i'^vi branching fractions using the isospin 
relation V{B^ T^-l+ve) = 2V{B+ TT°(.+vt) and the 



15|. The signifi- 



lifetime ratio tb+ /tbo = 1.071 ± 0.009 
cance of the B'^ — » r]t^vi signal is 3.2cr. 

A Bayesian 90% confidence limit B{B^ —* rj'l^vi) < 
OAT X 10^** is determined, assuming a fiat prior in the 
physical (nonnegative branching fraction) region, via the 
integral of the likelihood function from the signal extrac- 
tion, smeared by a Gaussian resolution function with 
varying width representing all other sources of uncer- 
tainty. We also determine the partial branching fraction 
AB{B+ ri'£+Vi) < 0.37 x 10"^ for < IGGcVVc^ 
and the ratio B{B+ ^ r]' i+vi) / B{B+ ^ 'qi+vi) < 0.57 
with 90% confidence level, the latter of particular im- 



portance in constraining the dynamics of the 77-77' system 
6]. These are in disagreement with a recently published 
resuh [H]. 

Extraction of \Vub\ from the measured B —> irivi 
branching fractions AB proceeds through the relation 
\Vub\ = ^/AB/{tboAC), with tbo = 1.530 ± 0.009 ps-i 
the meson lifetime 15] and the calculated reduced 
(i.e., appropriately normalized) decay rate over the cor- 
responding q^ range, which depends on the decay form 
factor Several form factor calculations are available, 
including one using light-cone sum rules ^121] and various 
lattice QCD methods 
Table im 
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Results are given in 
The branching fractions B{B r]^'Hi'£) will 
provide additional means of determining \Vub\ as accu- 
rate calculations of become available. 

In conclusion, we have measured the branching frac- 
tions for B — > Pive, where P are charmless pseudoscalar 
mesons, as a function of the squared momentum transfer 
q^. We report the total branching fractions, the third 
with a significance of 3.2a: 

B{B^ -> TT^i+i^i) = (1.38 ± 0.21 ± 0.07) x 10-^ (3) 
B{B+ TT°i+i^e) = (0.96 ± 0.15 ± 0.07) x 10"^ (4) 
B{B+ r]e+i^i) = (0.64 ± 0.20 ± 0.30) x 10"^ (5) 

with the first uncertainty statistical and the second sys- 
tematic, and, to 90% confidence level. 



B{B^ 



ri't+vt) < 0.47 X 10""*. 



(6) 



We combine the pionic branching fractions to obtain 



B{B" 



(1.54 ± 0.17 ± 0.09) X 10-"*, among 



the most precise measurements of this branching fraction 
available. We use the partial branching fractions to ex- 
tract \Vub\, using a variety of form factor calculations, 
and obtain values ranging from 3.6 x 10~^ to 4.1 x 10~^. 
The pionic branching fraction measurements represent a 
roughly 30% improvement over a previous BABAR mea- 
surement in this channel [s^, and is statistically indepen- 
dent of similar measurements in other channels 0, 0] . 
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FOM (The Netherlands), NFR (Norway), MES (Russia), 
MEG (Spain), and STFC (United Kingdom). Individuals 
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TABLE I; Partial and total branching fractions, in units of 10~ , for each decay channel; the first uncertainty given is statistical, 
the second is systematic. Ranges for are given in GeV^/c^. In the bottom row is the result from combining B*^ ~» -n' I'^u 
and B — > branching fractions. 

< 8 8 < < 16 > 16 < 16 total 

B" TT-i+u 0.59 ± 0.12 ± 0.03 0.34 ± 0.11 ± 0.02 0.46 ± 0.14 ± 0.03 0.92 ± 0.16 ± 0.05 1.38 ± 0.21 ± 0.07 

B+ ■k'^I+u 0.43 ± 0.09 ± 0.02 0.29 ± 0.08 ± 0.03 0.24 ± 0.09 ± 0.03 0.73 ± 0.12 ± 0.05 0.96 ± 0.15 ± 0.07 

B+ ril+v 0.28 ± 0.10 ± 0.01 0.16 ± 0.11 ± 0.01 0.21 ± Q.ntlxl 0-43 ± 0.15 ± 0.02 0.64 ± 0.20 ± 0.03 

B+ -> ri'i+v - - -0.05 ± 0.221 ji;!;^ 0.04 ± 0.22tg;g^ 

B" -> n'l+u (combined) 0.67 ± 0.10 ± 0.03 0.43 ± 0.09 ± 0.03 0.46 ± 0.11 ± 0.04 1.08 ± 0.13li');'og 1.54 ± 0.17 ± 0.09 



TABLE II: Values of \Vub\ derived using branching fractions 
measured in this Letter and various form factor calculations. 
Range for g^ is stated in GeV^/c^, reduced decay rate in ps~^. 
The given uncertainties on \Vub\ are, respectively, statistical, 
systematic and due to uncertainties in form factor calculation. 









AC 




(10-^) 


Ball & Zwicky \12] 


< 


16 


5.44 ± 1.43 


3.6 ±0.2 




Gulez et al. [17] 


> 


16 


2.07 ±0.57 


3.8 ±0.4 


±0.2t-I 


Okamoto et al. \18} 


> 


16 


1.83 ±0.50 


4.0 ±0.5 


±0.2«-i 


Abada et al. [19] 


> 


16 


1.80 ±0.86 


4.1 ± 0.5 


± 0.21 J:? 
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